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Abstract: We present updates of our results for neutral B-meson mixing and leptonic decay constants 
obtained in the quenched approximation from a mean-field-improved Shcikholeslami-Wohlert action 
at two values of lattice spacing. We consider quantities such as Bs d .. , /d (s) , fB (s) and the full AB = 2 
matrix-elements, as well as the corresponding S'f7(3)-breaking ratios. 



1. Introduction 

The study of Bd — Bd oscillations enables mea- 
surement of the poorly known CKM matrix ele- 
ment |Vtd|- The frequency of these oscillations is 
determined by 



Am d 



Ml 



Mi 



1.1 



where Mg and M^ are the heavy and light 
mass eigenvalues of the mixing system. Am, is 
experimentally measurable from tagged Bd me- 
son samples, and is also calculable in the Stan- 
dard Model. To leading order in l/Myp, the 
Standard Model prediction for Am, is |], ^ 



Am d 



^M^\V t dV t l\ 2 VB S (x t )C B (p) 



1 



2M B , 

.2 /7\,f2 



B d \0$ B = 2 (ri\B d )\ 



(1.2) 



where Xt — mf/Myy, i)Bi So(xt) and Cb(h) are 
perturbatively-calculated short-distance quanti- 
ties, ii is the renormalisation scale and Q^ B = 2 is 
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the four-quark operator [67'* (1— j 5 )d] [67^(1 — 
7 5 )d] . Since \V t b\ is equal to unity to very good 
accuracy, a measurement of Amd enables the de- 
termination of |Vtd|. The accuracy of this de- 
termination is currently limited by the theoret- 
ical uncertainty in the calculation of the non- 
perturbative strong-interaction effects in the ma- 
trix-element (Bd\ 0^ B=2 \Bd). An alternative 
approach [3j, in which many theoretical uncer- 
tainties cancel, is to consider the ratio, Am s / Amd, 
where Am s is the mass difference in the neutral 
B s — B s system. In the Standard Model, one has 



Am s 
Am d 



V td 
V ts 



Id 



M B . 
M Bd 

M Bd 
M B , 



Vtd 



M Bd 



M B . 



r sd 



(BJO 



AB=2 



{B d \Of B - 



\B d ) 



,(1.3) 



where Of ,B=2 is the same operator as £)f B=2 
with d replaced by s and where we have omitted 
the renormalisation-scale dependence of these op- 
erators as it cancels in the ratio. Because the uni- 
tarity of the CKM matrix implies |V( S |~|T4i,| and 
because \V c b\ can be accurately obtained from 
semileptonic B to charm decays, a measurement 
of Am s /Am d determines \Vtd\- This is a chal- 
lenging measurement and, at present, only a lower 
bound on Am s /Am d exists dJ. 

The matrix elements which appear in Eq. 
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( |l.3| ) are traditionally parameterised by 

MbM = (B q \0$ B = 2 (p)\B q ) 

= \M%J%B B ^) , (1.4) 

where q = d or s, where the -B-parameter, B Bq , 
measures deviations from vacuum saturation, cor- 
responding to B Bq — 1, and f B is the leptonic 
decay constant: 

(0167^7591^^)) = ip fi f Bq ■ (1-5) 

With this parameterisation, the quantity £ de- 



fined in Eq. (1.3) is given by 




fs d V B Bd 



(1.6) 



Because M Bs and M Bd are measured experimen- 
tally, £ 2 is the quantity in Eq. (1.3) which re- 



quires a non-perturbative determination. 

We report on B Bd(s) , f B(s) , B B J B Bd , f B J 
f B , r s d and £. Results for D-meson decay con- 
stants /d w and the SU(3) breaking ratio f B J 
ju are also given. These results are updates of 
those we presented in ||, |(| . 

2. Main features of the lattice calcu- 
lation 

Numerical calculations are performed in the quen- 
ched approximation at two values of the coupling, 
(3 = 6.2 and f3 — 6.0, corresponding to an inverse 
lattice spacing 1/a ~ 2.5 GeV (finer) and l/a ~ 
2.0 GeV (coarser), respectively. We use a mean- 
ficld-improved Shcikholcslami-Wohlcrt (SW) ac- 
tion |?| to describe the quarks. With this action, 
discretisation errors are formally reduced from 
0{a) to 0{a s a) and may be numerically smaller 
because of the mean- field improvement. This re- 
duction of discretisation errors is important in 
lattice calculations involving heavy quarks, be- 
cause these quarks have small Compton wave- 
lengths. For details on the parameters used in 
the numerical calculations, please refer to table 
H in Appendix [a|. 

At each value of the lattice spacing, we have 
three light quarks with masses in a range between 
~ m s /2 and ~ m s , which allows us to linearly 



extrapolate the quantities we are after to van- 
ishing quark mass and interpolate them to the 
strange-quark mass. To obtain results for the b 
quark, while keeping discretisation errors under 
control, we work with five heavy-quark masses 
straddling the charm mass 1 and extrapolate up 
to the b mass, guided by HQET. 

We use two methods to calculate r s( i: 

• Direct method: r s d is obtained from the 
direct calculations of M Bs and M. Bd - 

• Indirect method: r s d is obtained by cal- 
culating f B Jf Bd and B B jB Bd , and then 
combining them with the experimental M Bg / 
M Bd . 

We find that both heavy-quark-mass and light- 
quark-mass extrapolations are under better con- 
trol for the indirect method than they are for the 
direct method. 



3. Matching to the continuum and 
running in the MS scheme 

Results of lattice-regularised calculations have to 
be matched to the continuum renormalisation 
scheme in which Wilson coefficients are calcu- 
lated. We perform this matching at one loop 
||, ^, [l0|, [ll[] with mean-field improvement |l2| . 
At this order, it is consistent to use the tree- 
level value for improvement coefficient csw (see 
Appendix [X]). This is the procedure we use to 
obtain the central values for our results. 

Moreover, because chiral symmetry is explic- 
itly broken by Wilson fermions, the axial vector 
current requires a (multiplicative) renormal- 
isation and is related to its continuum counter- 
part via 

A™ nt = Z A {a s )Af\ (3.1) 

where Za is finite. 

For the four-quark operators, to subtract the 
contributions arising from the explicit breaking 



1 However, only three of these are used on the finer 
lattice (/3 = 6.2) when calculating matrix elements and 
_B-parameters of the four-quark operators. 
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of chiral symmetry, it is sufficient to consider the 
following basis of parity-conserving operators 

01,2 = 7 M x 7 M ± 7^75 x 7^75, 
03,4 = /x/± 75 x 75 , (3.2) 
O5 = u^v x a^y, 

where we only show their Dirac structure for sim- 
plicity. 0\ is the parity-even part of 0^ B=2 . 
This operator, in the MS scheme at the renor- 
malisation scale fi, is related to the above lattice 
operators by 

Or(^)-^nK,«M)Oi att (a) , (3-3) 

where 

5 

Oi att (a) = 0[ aU (a) + Y,Zu(a s )Of u (a) . (3.4) 

i=2 

Z\\ depends logarithmically on a/i. The Zu, i ^ 
1, account for the operator mixing due to the 
explicit chiral symmetry, and do not depend on 
a/i. 

The scheme of a s is not fixed at one loop. We 
choose a s = a-^ obtained from the procedure de- 
scribed in fll2| , with rif = (quenched approxi- 
mation) , which was shown to lead to particularly 
convergent perturbative expansions jL2| . Central 
values are obtained by identifying the scale of the 
coupling with the matching scale \i and match- 
ing at n = 2/a-a typical lattice ultraviolet scale. 
Running in the MS scheme is performed at two 
loops with the same coupling constant as for the 
matching, and rif = 0. 

4. Scaling with heavy-quark mass 

To study the behaviour of the various quantities 
with heavy-quark mass, we define 

$ KMp)^V^(^) 2//30 (4.D 
Z A \a s (M B ) J 

<M/i, M P ) = Bpfa) f^Mfl) 2 ^ (4.2) 
\a s {M B ) J 

(4.3) 



where Mp is the heavy- meson mass and Aip is 
the AF = 2 matrix element calculated at that 
mass. (3q is the one-loop coefficient of the QCD 
(3- function, with n./ = 0. In <&/, <3>b and $af=2, 
we have cancelled the logarithmic dependence of 
fp, Bp and Aip on Mp at leading- log order |13"| . 

For X(M P ) = <5> f , $ B , $af=2 and the SU{3)- 
breaking ratios, we use the HQET-inspired rela- 
tion, 

(4.4) 

to investigate the heavy-quark-mass scaling be- 
haviour of these quantities, as shown in figures 
0| and| 

5. Systematic uncertainties 

Our main results at the two values of lattice spac- 
ing are summarised in table [l]. In this table, the 
first error bar for each quantity is statistical. The 
other errors are systematic and we discuss them 
now. 

5.1 Discretisation errors 

In table |[ results for the decay constants dis- 
play significant variation with lattice spacing 2 . 
This suggests that discretisation errors for these 
quantities may be important. To quantify these 
errors we estimate residual, 0(arriQa s ) discreti- 
sation effects, associated with the mass rriQ of 
the heavy quark, as described in Appendix [b|. 
This is the second error bar on the decays con- 
stants and their S'L/(3)-breaking ratios, fn 3 / Jd 
and JeJIb- 

Such an estimate could, in principle, be car- 
ried out for the B-parameters and the SU(3)- 
breaking ratios r S( j and £. However, many of 

2 This poor scaling is not fully understood. It could be 
improved by using f n instead of m p to set the scale (see 
table |). For instance, f B at (3 = 6.0 would be ~ 191 MeV 
instead of 201 MeV while its value at f) = 6.2 would be 
160 MeV. However, m p is a valid means of setting the 
scale in quenched calculations and, as discussed below, 
we include a systematic associated with the uncertainty in 
the lattice spacing. Furthermore, this poor scaling is not 
present in the B-parameters and S{7(3)-breaking ratios 
which are the main thrust of the present work. 
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Figure 1: Scaling of <J?/ with heavy-quark mass on (a): the finer lattice and (b): the coarser lattice. The 
points labelled Mq s and Mqj, correspond to the heavy quarks, Q, used in our simulation. The curves are fits 
to the RHS of Eq. (4.4 ). The other points are the result of interpolation to Q = c or extrapolation to Q = b. 



these discretisation effects cancel trivially in the 
ratios of matrix elements defining these quanti- 
ties. Furthermore, a full quantification of 0(a 
rriQa s ) effects for Ms q and their U-parameters 
would require one to consider the mixing of the 



four-quark operators in Eq. (3.3) with operators 
of dimension seven, which is beyond the scope of 
the present work. Finally, in table [I], results for 
£>-parameters and their S'[/(3)-breaking ratios 
exhibit very little lattice-spacing dependence, sup- 
porting the idea that discretisation errors for these 
quantities are small. Thus, we assume that the 
statistical error for these quantities encompasses 
possible residual discretisation errors. For r ™ dlrect 
and £ , however, which are obtained using fs s / Jb , 
we take into account the discretisation error on 
this quantity. 

5.2 Matching uncertainties 

To estimate the systematic errors arising from 
the perturbative matching in the £>-parameters, 
we match at different \i in the range between 1 /a 
and it /a 3 , then run the resultant ^-parameters 
to 2/ a to compare them with the ones matched 
"directly" at 2/ a. The range [l/a,n/a] covers 
typical lattice ultraviolet scales and is vindicated 
by our study of Bk [H, where we find that 



continuum chiral behaviour is restored for these 
scales. We also consider the variation coming 
from computing Z\\{jjl = 2/ a) and Zu with the 
constant csw set to its mean-field-improved value 
instead of 1. All of these variations, which affect 
Bb and Bb 3 significantly, but not Bb s / Bb, are 
reflected in the second error bars on these B- 
parameters. 

Decay constants are independent of renor- 
malisation scale. However, the above procedure 
results in a ~ 4% change in Za through the 
/^-dependence of aiss(^t) and the change in the 
value of csw- This is reflected in the decay con- 
stants' third error bar but does not affect the 
corresponding 5C/(3)-breaking ratios. 

5.3 Heavy-quark-mass extrapolations 

As shown in figure [j], the decay constants have a 
pronounced extrapolation in heavy-quark mass, 
and the term quadratic in 1/Mp on the RHS of 



3 We always identify the scale of o^jg- with the matching 
scale. 



Eq. ([hj) contributes significantly. To quantify 
the systematic error associated with this extra- 
polation-the fourth error bar on the decay con- 
stants-we perform a fit of the heaviest three points 
in figure [I] to the RHS of Eq. ((P)), without the 
quadratic term. 

Figure ^ indicates that the linear heavy-quark- 
mass extrapolation of the i3-parameters works 
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Figure 2: Scaling of $s with heavy-quark mass on (a): the finer lattice and (b): the coarser lattice. The 
points labelled Mq s and Mqa correspond to the heavy quarks, Q, used in our simulation. The curves are fits 
to the RHS of Eq. (4.4) without the quadratic term in 1/Mp. The other points are the result of extrapolation 
toQ = b. 



well and is mild: the associated uncertainty should 
be covered by the statistical error. 

The AF — 2 matrix elements have a very 
pronounced dependence on heavy-quark mass, as 
seen in figure |^. Since we are not reporting re- 
sults for the individual MB d{s) , we do not quan- 
tify the systematic errors associated with their 
determination. However, this strong mass-de- 
pendence is one of the elements which make a 
reliable determination of r sc i, from the ratio of 
individually calculated M.B d{s v difficult p5| . 

5.4 Uncertainties in the determination of 
the lattice spacing 

In quenched calculations, the value of the lat- 
tice spacing varies significantly with the quantity 
used to set the scale. This variation may be due 
to quenching effects, as well as any other system- 
atic uncertainty which may affect the quantity 
used to set the scale. In this work, we deter- 
mine the lattice spacing from the p-meson mass 4 . 
We then vary the inverse lattice spacing, 1/a, in- 
creasing it by 10% and decreasing it by 5%. This 
range covers the typical variations observed in 
the determination of the scale from gluonic or 



light-hadron spectral quantities, for the action 
and parameters we use |^ . 

Uncertainties in the lattice spacing will ob- 
viously affect the determination of all the de- 
cay constants, as they are dimensionful. They 
will also slightly change the curves in the heavy- 
quark-mass extrapolations (figures [l], || and ||). 
Furthermore, they induce a variation in the stran- 
ge-quark mass, which we obtain from the mass of 
the kaon, and therefore affect all quantities which 
depend on this mass. 

In practice, we find that the variation of the 
lattice spacing discussed above does not induce 
a significant change in the B-parameters. How- 
ever, it does affect all the decay constants and 
5'J7(3)-breaking ratios. This is reflected in the 
last error bar on these quantities. 

5.5 Quenching errors 

Quenching errors have been studied using quen- 
ched Chiral Perturbation Theory (q^PT) and have 



been found to be small for the B-parameters [17 



4 The scale determined from gives compatible re- 
sults at (3 = 6.2, as shown in table ^. 



Moreover, numerical simulation with two light 
flavours of dynamical quarks indicate that those 
in fs s I Jb &re also small |J, |l). Thus, they 
ought to be small for r sc t and £. Contrary to 
this, quenching errors may be significant for the 
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Figure 3: Scaling of $af=2 with heavy-quark mass on (a): the finer lattice and (b): the coarser lattice. The 
points labelled Mq s and Mq^ correspond to the heavy quarks, Q, used in our simulation. The curves are fits 
to the RHS of Eq. (4.4 ) (without the quadratic term in 1/Mp for the finer lattice). The other points are the 
result of extrapolation to Q = b. 



decay constants themselves, as indicated by both 
q%PT |l7j and numerical simulation |D| Q . 

As we mentioned in the previous section, the 
uncertainty in the lattice scale may be, in part, 
a quenching effect. Thus, to the extent that it 
is, we have already accounted for some quench- 
ing errors. A more thorough estimate of these 
effects, however, would require a dedicated un- 
quenched simulation which is beyond the scope 
of this work. Therefore, we do not attempt to 
quantify quenching errors any further. 

6. Final results 

Since two lattice spacings are not sufficient for 
an extrapolation to the continuum limit (a = 0), 
we quote the results obtained on the finer lattice 
([3 = 6.2) as our best estimates. And because it 



appears to be more reliable, we quote 



^indirect 



for 



r s d- 



Our main preliminary results are thus 



£ = 1.15(6)+ 2 



3 ) 



r sd = 1.37(14)12 , 



= 0.98(3) , 
= 1-16(6)1 2 3 , 



B B (5 GeV) = 0.92(4) lg , 
B Bs {5 GeV) = 0.91(2)1.3 , 
!b = 161(16)1?| MeV, 
f B . = 192(14)1 2 ! MeV , 
f D = 195(10)1 22 MeV, 



-10 
1-21 



MeV 



= 1-15(4)11 , 



f D , = 224(7)"_ 9 
/£» 

fD 

where the first error bar is statistical and the 
second is systematic, the result of adding our 
long list of systematic errors in quadrature. By 
choosing the results obtained on the finer lattice, 
we are also being conservative in our estimate of 
statistical errors as they are larger than on the 
coarser lattice for which we have higher statistics. 

Note that our results, obtained with large 
statistics and a highly improved action are com- 
patible with recent world averages || |l^, ||l) . 
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lattice 


coarser 


finer 


3 


6.0 


6.2 


fn (MeV) 


^ UJ -V°/— 0— 4—0- 8 


224(7)+"+*+"+^ 

^^ny i y_o-3-0— 8 


fn(MeV) 


&A'±y±j_ _ 4-0-10 


195(10)+£+I+^ 

± jo\±\j) _o-3-0-10 
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f D Jf D 
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Ib 3 1 Jb 


i.i4(2)^i i + :i 


1.16(6)^1 i i 2 3 


Bb 3 (5 GeV) 


0.92(2)+ +4+ +° 


0.91(2)1 t 3 t ±° 


B B {5 GeV) 


0.90(4)1 t 4 t +° 


0.92(4)+ + 3 + +° 


B B jB Bi 


1.02(3)1 t° t t° 


0.98(3)+ +° + +° 


^indirect 
' sd 


1.38(7)+^+ 1 1 4 6 


1.37(14)^1 1 + J 6 


direct 
' sd 


1.52(18)+ + + +% 


1.71(28)+ + + + £ 




1.15(3)±£ t ?i 


1.15(6)1 1 1 1- 



Table 1: Results at the two values of lattice spacing. rf d dilcct = (t^^) 2 "^ and r^ rcct = (A4 Ss /.MB d )- 
The first error bar on each quantity is statistical while the others are systematic, as described in Section ^. 
Blank error bars are put in to help keep track of which systematic effect each error corresponds to. 
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A. Simulation Details 

The quarks in the simulation are described by 
the SW action 



S sw 



(A.l) 

where S^p is the standard Wilson action, go the 
bare gauge coupling, a lattice realisation of 
the Yang-Mills field strength tensor, n q the ap- 
propriate quark hopping parameter and csw ; the 
so-called clover coefficient, csw = 1 corresponds 
to tree-level improvement. With this value of 
csW) leading discretisation errors are 0(a B a) in- 
stead of 0(a) as they are with the standard Wil- 
son action, corresponding to csw = 0. We actu- 
ally use a mean-field-improved SW action with 
values of csw given in table g, where the param- 
eters used in our simulations arc summarised. 



lattice 


coarser 


finer 




6.0 


6.2 


lattice size 


16 3 x 48 


24 3 x 48 


csw 


1.47852 


1.44239 


# of cfs. 


498 


188 


a-^Mp) (GeV) 


1.96(5) 


2.54(8) 


a- x {U) (GeV) 


1.87(4) 


2.52(8) 



Table 2: Simulation parameters. a^ 1 (M p ) and 
a~ (/ir) &re the values of the inverse lattice spacing 
determined from calculations of the p-meson mass 
and the pion decay constant, respectively. The lat- 
ter is given for u, — 2/ a. 



B. (9(a)-improvement of the axial cur- 
rent 

The leading discretisation errors with the mean- 
field- improved SW action are formally of 0{a s a) , 
as they are for the tree-level improved SW ac- 
tion. To estimate these errors, we consider the 
following variation in our procedure. 

C(a s a)-improvement of the axial current re- 
quires one to include the effect of the ad^P (P 
the pseudoscalar density) counterterm through 
the replacement 



An + c A aduP , 



(B.l) 
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as well as to rescale the quark fields as 



7] 7 h A — 

(1 + —am q )q 



(B.2) 



with both ca and 6^ evaluated at one loop ||22|~ 
\23A. Thus, from a comparison of results obtained 



[9 



with ca and set to their tree-level values (ca = 



(J and bA — 1) to those obtained with ca and 
6a evaluated at one loop, we can have an esti- 
mate of the effect of 0(a s a) discretisation errors. 
We do not use the one-loop results as central 
values for the decay constants, for consistency 
with nnr determination of the ^-parameters. Tn- 



[10 

[11 

[12; 



rleerl; (r7(rv J a)-improvement, nf the fnur-gnark op- 



erators would require one to consider the mixing; 



[13 



of these operators with operators of dimension 



seven, which is beyond the scope of the present [14 
work. 



To correct for some higher-order discreiisa- 

tion effects, we actually use KLM normalisation 
PH for the quark fields. Thus, our central values 
are obtained with the normalisation 



7/1 + am q q 



TTT3T 



[15; 

[ie; 

[17 



and tho nno-lnnp variarir 



[18 



1 + arriq/2 



(1 



, l-loop 
'A 



-am q )q . (B.4) 



We also check that results obtained with tree- 
level normalisation Q)a 
within the diacrctiaation error we quote. 



1 in Eq. (B.2)) p§ lie 



[19 



[20 
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